Role of External Flow and Frame Invariance in Stochastic Thermodynamics 
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For configurational changes of soft matter systems affected or caused by external hydrodynamic 
flow, we identify applied work, exchanged heat, and entropy change on the level of a single trajectory. 
These expressions guarantee invariance of stochastic thermodynamics under a change of frame of 
reference. As criterion for equilibrium vs. nonequilibrium, zero vs. nonzero applied work replaces 
detailed balance vs. nonvanishing currents, since both latter criteria are shown to depend on the 
frame of reference. Our results are illustrated quantitatively by calculating the large deviation 
function for the entropy production of a dumbbell in shear flow. 
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Thermodynamic notions like applied work, dissipated 
heat, and entropy have been used successfully to an- 
alyze processes in which single colloidal particles or 
biomolecules are manipulated externally Various ex- 
act relations have been shown to constrain the distribu- 
tion function arisingfrom the ever present thermal fluctu- 
ations on this scale [2I, [^, 3, 01 ■ So far, the systematic con- 
ceptual work has been focused on cases where the driving 
crucial to generate a nonequilibrium situation arises from 
a time-dependent potential expressing the effect of a mov- 
ing laser trap, micropipet, or AFM tip. In these cases, 
the identification of external work, internal energy and 
hence dissipated heat is straightforward. Many experi- 
ments on a variety of systems have proven the consistency 
and potency of such an extension of thermodynamic no- 
tions to the micro or nano world [1, .7, 8, 9, fol. 

As another source of nonequilibrium, external hydro- 
dynamic flow arguably is the most common and best 
studied case in soft matter systems. On the level of sin- 
gle objects like a polymer, vesicle, or red-blood cell, it 
can cause dramatic shape transitions (for recent exam- 
ples see, e.g., Refs. [H, [H, Ei, El, [ll). The theoret- 



ical analysis of these phenomena is typically based on 
equations of motion like the Langevin equation for sin- 
gle objects or distribution functions and their projection 
to few-body correlation functions for bulk systems like 
colloidal suspensions Concepts like work or 

entropy production, however, have played no systematic 
role in describing such phenomena yet. Having in mind 
the conceptual advantage achieved for time-dependent 
potentials by using such notions, the question arises quite 
naturally whether a similar analysis for nonequilibrium 
phenomena in soft matter systems caused by external 
flow is possible. Indeed, Turitsyn et al. [lO] calculated the 
entropy production for linear equations of motion and il- 
lustrated it for a dumbbell in shear flow. However, their 
identification of entropy production was not related to 
the dissipated heat leading to a fundamentally different 
dependence on the external flow. 

Conceptually, there is an even deeper issue hidden be- 
hind the proper treatment of external flow. Its presence 
or absence depends on the frame of reference as the pres- 



ence or absence of currents does. Since the latter are usu- 
ally taken as an indicator for nonequilibrium, our analysis 
will question the traditional role of both detailed balance 
and nonvanishing currents as criteria for equilibrium and 
nonequilibrium, respectively. They will be replaced by a 
frame-independent identification of applied work. 

For an almost trivial but revealing paradigmatic case, 
consider a colloidal particle dragged through a viscous 
fluid along the trajectory x{t) by a harmonic laser trap 
of strength k moving with velocity uq [^, @| . In the labo- 
ratory frame, the particle is moving in a time-dependent 
potential U{x,t) = {k/2){x — uot)^. By applying the 
standard definition of work Q as the external change of 
the potential energy, the applied power 



^dtU 



(1) 



becomes w — ~~UQk{x — uqI). Changing to the comoving 
frame with y = x — u^t, the potential U{y) — (fc/2)?/^ be- 
comes time-independent and therefore one would naively 
find w = 0. Of course, the work should be independent 
of the chosen frame of reference. Moreover, since in the 
comoving frame detailed balance holds and no current 
in the y variable occurs, one would usually consider the 
system to be in equilibrium in this representation. The 
physical origin of the apparent contradictions is the fact 
that in the comoving frame the particle experiences a 
steady flow of the fluid with velocity —uq. Due to fric- 
tion, this flow pushes the particle against the force of the 
potential thus spending work. The standard definition of 
work fails in the presence of flow since the flow advects 
the particle, which is not accounted for in Eq. ([T]). Modi- 
fying this expression will affect also the expression for the 
heat dissipated into the fluid as well as the expression for 
the entropy production. 

We will now derive the refined expressions for work 
and heat in the presence of a, possibly time-dependent, 
external flow u(r, t). For complete generality and future 
applications, we consider a soft matter system like poly- 
mers, membranes, or a colloidal suspension composed of 
N particles with positions F = {ri, . . . , tat}. The index 
denotes the particle number and in the following we sum 
over same indices. The system has a total energy U{r, t) 



2 



which is the sum of an internal energy due to particle 
interactions and a possible contribution due to external 
potentials. In addition to direct interactions, we allow for 
hydrodynamic interactions which are of great importance 
for soft matter systems. 

In 1997, Sekimoto formulated the first law of thermo- 
dynamics 6q = Sw — dU along a single stochastic tra- 
jectory [2lj. We follow this route and start by identi- 
fying three possible sources of work w, i.e., of changes 
of energy caused externally. First, the potential energy 
C/(r, t) can be time-dependent if it is manipulated exter- 
nally. Second, it is convenient to allow for direct noncon- 
servative forces ft that spend work through displacing 
the particles. Third, the external flow drives the sys- 
tem and therefore changes its energy, which must be 
taken into account as work as well. The advection of 
particles is described through the convective derivative 
Dt = dt + u(rfc) • Vfc. Through both replacing the par- 
tial derivative in Eq. ([1]) with the convective derivative 
and measuring particle displacements with respect to the 
flow, we obtain the new definition 



w 



BtU + tk- [i-fe -u(rfe)] 



(2) 



for the work increment 6w = ibdt which replaces Eq. ([T]) . 
The first law then leads to the heat production rate 

q = w- dU/dt = [Vk - u(rfe)] ■ [-VfeC/ + fk]. (3) 

The sign of heat is convention, here we take it to be 
positive if energy is dissipated into the surrounding fluid. 
The expression ([3]) shows that no heat is dissipated when 
particles move along a fluid trajectory, = u(rfe). In 
this case, the applied work corresponds to the change in 
total energy. 

A change of frame necessarily changes the flow there- 
fore leaving the expressions for work ^ and heat ^ in- 
variant. Explicitly, we can allow time-dependent orthog- 
onal transformations and an arbitrary time-dependent 
shift of the origin. Coming back to the introductory ex- 
ample, we see that for a stationary trap in the comoving 
frame with flow u = — mq, the applied power following 
Eq. ([2]) now reads w — —u^ky in agreement with the 
rate obtained for a moving trap and resting fluid. Finally, 
note that the definitions ^ and ^ are independent of 
whether or not hydrodynamic interactions arc induced. 
These will affect the dynamics but they do not enter the 
definitions of work and heat explicitly. 

We now turn to entropy production. A trajec- 
tory dependent stochastic entropy is defined as s{t) = 
— ln^/;(r(t), i) through the distribution function ijjiT^t) 
of the N particle positions The temperature of the 
surrounding fluid is T and throughout the paper, we 
set Boltzmann's constant to unity. Calculating the to- 
tal time derivative s = ds/dt, we obtain the equation of 



motion 

s = dts-vk- Vfc In^ = Dts + [vk - u(rfc)] • Tk/T 

-[rfc-u(rfc)]-[-VfcC/ + ffc]/r, (4) 

where we have separated the heat production rate ([3]) 
and introduced the total effective force 



-Vk{U + T\ni,)+h 



(5) 



Apart from the gradient of the potential energy J7(r,t) 
and nonconservative forces ffc that cannot be written as 
gradient of a potential, a "thermodynamic" force con- 
tributes to J-k arising from the stochastic interactions 
between system and the surrounding fluid 



17|. 



In the 

absence of external flows and nonconservative forces, de- 
tailed balance must hold. The thermodynamic force 
then ensures that Vk{U + Tin?/') = leads to the cor- 
rect Gibbs-Boltzmann equilibrium distribution ipcq ~ 
exp(-C//r). In Eq. (g]), the last term is the heat © 
divided by the temperature of the fluid T. We interpret 
this term as the entropy produced in the fluid through 
Clausius' formula Asm = l/T. The total entropy pro- 
duction rate stot = s + Sm then becomes 



Stot 



Bts+[rk-u{rk)]-Tk/T. 



The mean total entropy production rate follows as 
T(stot) = ([rfc -u(rfc)] •:Pfc) 



(6) 



(7) 



since the mean of the convective derivative is zero for 
incompressible fluids (V • u = 0) and vanishing boundary 
terms. 

So far we did not resort to a specific dynamics. How- 
ever, in order to both prove positivity of Eq. ([7]) and give 
it a more familiar appearance known from the theory 
of polymer d yna mics [17| , we turn to the Smoluchowski 



equation jl7l. Il8l| 



5tV + Vfc • (vfcV) = 



(8) 



governing the evolution of the distribution function 
TpiT, t). Any deviation of a particle's local mean velocity 



Vfc = u(rfc) + n^^iTi 



(9) 



from the velocity of the external flow u(r), which drags 
the particles due to friction, has to be caused by the total 
force J-k ■ In this model, hydrodynamic interactions now 
enter through a dependence of the symmetric mobility 
matrices t^^ki^X) on the particle positions F. We demand 
the inverse matrices defined through = l5ki to 

exist. Specific expressions for /.tj.; are obtained from ei- 
ther the Oseen or Rotne-Prager tensor [l8l |. 

The local mean velocity Vfc(r,t) is the average of the 
actual stochastic velocity rfc over the subset of trajecto- 
ries passing through a given configuration F. This allows 
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us to perform the mean in Eq. ^ in two steps: First, 
we average over stochastic trajectories crossing a specific 
point r in configuration space, which amounts to replac- 
ing Tfc by Vfc. Second, we average over the distribution 
ip{T,t) of the point F. FinaUy, we use Eq. ^ to replace 
the total force J^k- The resulting quadratic expression 

r(itot) = (K - u(rfe)] ■ Mfc/[vi - u(rO]) ^ (10) 

is well known phenomenologically as the energy dissipa- 
tion function in the theory of polymer dynamics 17]. It 
is nonnegative with the equal sign holding in equilibrium 
only where JF^ = 0. 

On the trajectory level, we have identified heat and 
entropy production using physical arguments. There is a 
more formal way of identifying entropy by starting from 
the weight P[F(t)|Fo] ~ exp{-5[F(i)|Fo]/T} for a tra- 
jectory F(i) starting in F(0) = Fo involving the action 
S. The entropy produced in the surrounding medium 
is usually identified with the part of the action that is 
asymmetric with respect to time- reversal (2^ . For ex- 
plicit expressions for both the path probability and the 
action in the presence of hydrodynamic interactions as 
needed in our case, we refer to Ref. ,23]. 

We can deduce two expressions for the asymmetric 
part of the action since in the presence of external flow, 
we have two choices how to define the operation "time 
reversal". First, we could formally treat the flow as a 
nonconservative force. It would then be invariant with 
respect to time reversal and we would obtain 



ik ■ [-^kU + fk + f^kMri)]/T. 



(11) 



This expression has been identified with the entropy pro- 
duction rate in Ref. Second, we can extend the 
operation of time reversal to the particles of the fluid, 
which is physically more appropriate. Reversing the ve- 
locity of the fluid particles then effectively amounts to 
the change of the sign of the external flow velocity u(r). 
This leads to Sm = q/T involving the heat production 
rate from Eq. ([3]). We thus recover the physically moti- 
vated definition of the heat ^ as the asymmetric part of 
the action under time reversal only if the flow is reversed 
as wen (23. 

For a discussion of the crucial difference between 
Eqs. (fTT|) and ((3]) in a specific system, consider a dumb- 
bell in a flow u(r) = Kr with traceless matrix k = 

+ K^, which can be split up into a symmetric part 
and a skew-symmetric part k^. The dumbbell consists 
of two particles at positions ri and r2 with displacement 
r = Ti — r2 and potential energy U{r) ~ (A:/2)r^. The 
applied power ([2]) then becomes 



w 



u(r) • VC/ = fcr 



(12) 



i.e., only the symmetric part kP , which reflects an elonga- 
tion flow component, contributes to the work. Physically, 



the elongation flow drags the two particles apart due to 
friction and therefore spends permanently work against 
the elastic force keeping the particles together. Neglect- 
ing the boundary term AU/t '--^ in the long time limit, 
the time-extensive part of the entropy production rate 
is then Sm ^ {k/T)v ■ k^y. Note that here only the dis- 
placement r enters. This is not true if the expression pT|) 
is interpreted as entropy production rate, leading to the 
time extensive part 

sl,^{T^JiQ)~^[2Yl■K^Yl+{l/2)r■K\] (13) 

with bare mobility /iq. Here, also the center of mass 
R = (ri -|-r2)/2 contributes. This is obviously physically 
unsound as the center of mass undergoes free diffusion 
and therefore does not produce entropy in the medium 
on average. Moreover, even in the second term of Eq. (fT3|) 
concerning the relative motion, the source of dissipation 
is the skew-symmetric part in contrast to Eq. (jl2p . 

We will now calculate the large deviation function of 
the medium entropy production rate for the dumbbell in 
two-dimensional shear flow, i.e., all entries of the matrix 
K are now zero except Kxy = 7 where 7 is the strain 
rate 



181 ]. The large deviation function 



h{a) — lim — lnp(Asni,i) 

t — >oo t 



(14) 



quantifies the asymptotic fluctuations <j = As^/ {{sm)t) 
of the entropy production in the limit of large observa- 
tion times t with mean production rate (sm). Instead 
of obtaining the function h{a) from the time-dependent 
probability distribution p(Asni, t) directly, we will calcu- 
late its Legendre transform ao('^) as the lowest eigenvalue 
of the operator L\ — Lq ~ Aim- The operator L\ gov- 



erns the evolution of the generating function [2J, |25|, |26| . 
The production rate depends only on the displacement r, 
whose dynamics is determined through the Smoluchowski 
operator [13] = V^'^ ' + {T/k)V] - kv ■ V . For 
this illustration, we neglect hydrodynamic interactions 
by using mobility matrices /Xj,; = ^lo'^^m defining the 
time scale t^^ = 2/zofc. 

For the lowest eigenfunction of the eigenvalue equation 



L\ipo{v, A) = -ao(A)V'o(r, A), 



(15) 



we use the ansatz V'o = exp[— (fc/2r)r • Cxy] with a sym- 
metric matrix Ca- We are lead to this ansatz since we 
know that the stationary distribution (A = 0) of r is 
a Gaussian [23], where Cq becomes the inverse covari- 
ance matrix. Inserting this ansatz into Eq. (jl5[) results in 
ao(A) = Tq"^ tr(CA — 1), where Ca is the solution of the 
quadratic matrix equation 



(Ca + D)T(Ca+D) 

1 7(2A-1)- 
7(2A-1) 1+72 



^-2 



-1 7 
-1 
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FIG. 1: Left: Comparison of the eigenvalue ao{X) from 
Eq. dl) (solid) with the solution ^ obtained in Ref. HO] 
(dashed) based on Eq. (Parameters for both curves are 

7 = 3 and tq — 1.) The mean entropy production rate is the 
same for both expressions as indicated by the matching slope 
at A = 0. Right: The corresponding large deviation functions 
h{a) obtained from the Legendre transform of ao(A). 



The only parameter left is the dimensionless strain rate 
7 = 7T0. We only need the trace of the matrix Ca, which, 
after some tedious calculations, can be expressed as 



trCA = ytrSA + 2VdetSX-72/4-trD 
leading to the solution 

ao(A) = Tf^'^ 



^i + iyTT4^2^(l-A)-l 



(16) 



This function is defined within the interval A_ $J A ^ 
with branch points A± = i(l ± \/l + 7~^). Such branch 
points imply linear asymptotes for h{a) and therefore 
exponential tails for the distribution p(Asm) [31}. On the 
level of the function ao(A), the fluctuation theorem [H!, 
0j m, nil is expressed through the symmetry ao(A) — 
ao{l - A), which is fulfilled by Eq. In Fig. [H we 

compare the function p6p with the solution 



«S(A) 



v/l + 72A(l-A) -1 



(17) 



obtained in Ref. [20| based on the entropy production 
rate in Eq. with fixed center of mass (R = 0), which 
also fulfills a fiuctuation theorem. Note that our expres- 
sion for h{a) based on Eq. predicts substantially 
larger fluctuations especially for trajectories with larger 
than mean entropy production. 

The dumbbell just discussed also provides a counter- 
example to the standard definition of nonequilibrium. In 
pure rotational flow with skew-symmetric matrix 



K = H 



7/2 

-7/2 



there is a nonvanishing current manifested through the 
ever tumbling dumbbell superficially indicating nonequi- 
librium. On the other hand, the length |r| undergoes only 
equilibrium fluctuations. For all practical purposes, this 
is an equilibrium system. In fact, the applied power ^ 



vanishes. Combining this result with the moving trap 
case discussed at the beginning, we conclude that in the 
presence of flow the proper frame-invariant criterion for 
distinguishing equilibrium from genuine nonequilibrium 
is whether or not the applied power obeys w — 0. 

Summarizing further, we have shown that consider- 
ing external flow in the expression for the work ^ [and 
therefore in the derived quantities like dissipated heat ([3]) 
and entropy production ([6])] guarantees frame invariance 
of stochastic thermodynamics. While formally two ex- 
pressions for the asymmetric part of the path probability 
under time reversal are possible which both fulfill the 
fluctuation theorem, only one leads to a consistent iden- 
tification of thermodynamic quantities for the dynamics 
of soft matter systems in flow. 

We thank R. Finken and R. K. P. Zia for inspiring 
discussions. 
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